Dynamics of rumor propagation on small-world networks 
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We study the dynamics of an epidemic-like model for the spread of a rumor on a small-world 
network. It has been shown that this model exhibits a transition between regimes of localization 
and propagation at a finite value of the network randomness. Here, by numerical means, we per- 
form a quantitative characterization of the evolution in the two regimes. The variant of dynamic 
small worlds, where the quenched disorder of small- world networks is replaced by randomly chang- 
ing connections between individuals, is also analyzed in detail and compared with a mean-field 
approximation. 
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I. INTRODUCTION 

The networks that underlie real social interactions, 
whose nodes represent single individuals and whose links 
connect individuals that are expected to interact, vary 
with time and strongly depend on the kind of interactions 
involved. Generally, however, social networks exhibit two 
specific topological properties that are closely related to 
the nature of social interactions. First, they are highly 
clustered, which means that two randomly chosen neigh- 
bors of a given individual have a relatively large proba- 
bility of being in turn mutual neighbors. Second, the dis- 
tance between any two nodes in the network, measured 
as the number of links of the minimal path connecting 
the two nodes, is on the average very small as compared 
with the total number of nodes or links. This is the so- 
called small- world effect. Small- world networks (SWNs) 
constitute a mathematical model for social networks that 
captures these two properties They are partially 

disordered networks, interpolating between regular lat- 
tices and fully random graphs. In fact, iV-node regular 
lattices with connections beyond nearest neighbors have 
high clustering, but the average distance between nodes 
is of order N . On the other hand, the average distance 
in random networks is of order In A'' <g; A^, but the prob- 
ability that two neighbors of a given node are mutual 
neighbors is of order iV~^. For moderate disorder, SWNs 
preserve the two desirable properties of ordered and ran- 
dom networks (|-|^, and are therefore a convenient tool 
for the mathematical study of social processes . 

Small-world networks are built starting from an or- 
dered lattice with moderately high connectivity, which 
insures high clustering. Then, each link is removed with 
probability p and reconnected between two randomly se- 
lected nodes. This process creates a shortcut between 
two otherwise distant regions of the network. The prob- 
ability p measures the degree of disorder or randomness 
of the resulting graph. For p = Q order is fully preserved, 
while for p = 1 & random graph is obtained. Note how- 
ever that the average connectivity is constant. 



Topological properties of SWNs, such as the average 
distance between nodes, display a crossover from the be- 
havior observed in regular lattices to that of random 
graphs at a randomness p ^ N^^. In the limit of an in- 
finitely large network, a critical transition between both 
regimes occurs at pc = |@,||. Similar transitions at 
the same critical point are found for some simple dy- 
namical processes on SWNs, such as for Ising-like spin 
systems and ensembles of coupled oscillators. In con- 
trast, it has been recently shown that other kinds of pro- 
cesses display a transition between qualitatively different 
dynamical regimes at finite values of the randomness. 
Specifically, in an epidemiological model where an ini- 
tially susceptible individual infected by contagion under- 
goes a disease cycle that returns to the susceptible state, 
a transition at finite p occurs between a regime where 
the disease cycles of different individuals are temporally 
uncorrected (low p) to a regime where the cycles synchro- 
nize (highp) Moreover, in an epidemic- like model for 
rumor propagation a quantitatively similar transition has 
been found between a regime where the rumor remains 
localized (low p) to a regime where it spreads over a finite 
fraction of the network (high p) jl^ . 

The critical-phenomenon nature of the transition 
found in the model of rumor propagation has been con- 
vincingly proven by means of finite-size scaling analysis 
]lOt . This paper, on the other hand, focuses on a detailed 
characterization of the dynamical properties of the same 
model, with emphasis on the effects introduced by the 
small-world topology. In the next section we introduce 
the model and summarize the main results on the crit- 
ical transition between the regimes of localization and 
propagation. The core of the paper. Sects. Ill and IV, 
is devoted to establish the connection between the sev- 
eral parameters of our model and suitable quantities that 
characterize its evolution. This is done both in quenched 
small- world networks and in the so-called dynamic small 
worlds, where quenched disorder in the interaction links 
is replaced by stochastic choice of the interaction part- 
ners. We emphasize similarities and differences between 
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both cases. Finally, we summarize and discuss the main 
results. 



II. MODEL OF RUMOR PROPAGATION 

Consider a population formed by N individuals where, 
at each time step, each individual adopts one of three 
possible states. In the first state, the individual has not 
heard the rumor yet. In the second state, the individual is 
aware of the rumor and is willing to transmit it. Finally, 
in the third state, the individual has heard the rumor 
but has lost the interest in it, and does not transmit it. 
By analogy with SIR epidemiological models these 
three states are respectively called susceptible, infected, 
and refractory. At the beginning, only one individual is 
infected and all the remnant population is susceptible. 
The dynamical rules act as follows At each time 

step, an individual i is chosen at random from the in- 
fected population. This individual contacts one of her 
neighbors, say j. If j is in the susceptible state, i trans- 
mits the rumor and j becomes infected. If, on the other 
hand, j is already infected or refractory, then i loses her 
interest in the rumor and becomes refractory. 

In qualitative terms, the dynamics can be summarized 
as follows. In the first stage of the evolution, the number 
of infected individuals increases and, at a lower rate, the 
refractory population grows as well. As a consequence, 
the contacts of infected individuals between themselves 
and with refractory individuals become more frequent. 
Eventually, the infected population begins to decline and 
vanishes, and the evolution stops. At the end, the popu- 
lation is divided into a group of Nn refractory individu- 
als, who have been infected at some stage during the evo- 
lution, and a group of susceptible individuals who have 
never heard the rumor. It has been shown that, in 
the case where contacts can be established between any 
two individuals in the population, the fraction r = Nr/N 
of refractory individuals at the end of the evolution ap- 
proaches a well-defined limit r* for asymptotically large 
systems, N — > oo. This fraction is given by the nontrivial 
solution to the transcendental equation 

r* = 1 -exp(-2r*), (1) 

i.e. r* K, 0.796. In other words, some 20% of the popu- 
lation never becomes aware of the rumor. 

We are here interested in the case where contacts be- 
tween individuals are established along the links imposed 
by a socially plausible structure, namely a small-world 
network. As advanced in the Introduction, SWNs are 
built starting from an ordered lattice, with one individ- 
ual at each node. We choose a one-dimensional array 
with periodic boundary conditions, where each node is 
connected to its 2K nearest neighbors, i. e. to the K 
nearest neighbors clockwise and counterclockwise. Then 
each of the K clockwise connections of each node i is 
rewired with probability p to a randomly chosen node 



j, not belonging to the neighborhood of z. A shortcut 
is thus created. We avoid double and multiple links be- 
tween node, and discard realizations where the SWN be- 
comes disconnected. 

Previous analysis of this system, focused on the char- 
acterization of its final state, has revealed that a critical 
transition between two well-differentiated regimes occurs 
at a finite value Pc of the randomness For p < pc, 
the final number Nfi of refractory individuals, averaged 
over many realizations of the system, is independent of 
the population size. Therefore, as ^ oo, the ratio 
r — Nfi/N tends to vanish. In such limit, only an in- 
finitesimal fraction of the population becomes aware of 
the rumor, which remains localized in a small neighbor- 
hood of its origin. On the other hand, for p > pc the 
average value of r approaches a constant as N grows. 
Finite-size scaling analysis in the specific case of = 2 
shows that, for asymptotically large systems, 

r^\p-Pc\\ (2) 

with Pc ~ 0.19 and 7 « 2.2. For larger values of K, 
the critical randomness Pc decreases. The exponent 7, in 
contrast, seems to be universal. 

A clue to the origin of the localization-propagation 
transition is provided by the distribution /{Nn) of val- 
ues of Nn over large series of realizations of our system 
for fixed K, p, and N. Figure |l| shows those distribu- 
tions for K = 2, two values of p — below and above the 
transition — and three values of N. At each realization, 
the SWN is constructed anew. For p < pc the distribu- 
tion is approximately exponential, and does not depend 
on N. Consequently, the average value of Nfi^ is also in- 
dependent of the system size and, as advanced above, 
the ratio r = Nn/N decreases as TV grows. More specif- 
ically, r ^ N^^ for large N. In a typical realization for 
p < Pc the rumor remains localized due to the high in- 
terconnectivity of the network at the local level and the 
scarce density of shortcuts. Transmission occurs between 
a small group of individuals which rapidly loose their in- 
terest in the rumor, and propagation to distant regions 
is highly improbable. 

In contrast, the distribution for p > pc is bimodal, 
with two maxima separated by a local minimum. The 
small- Afl regime is still independent of N and attains a 
maximum near Nn = 0. On the other hand, for large 
values of Nn, we find an additional bump-like structure, 
which changes with the system size. Specifically, the po- 
sition of its maximum N'^'^^ shifts rightward as N grows, 
as N^^^ « 0.25A^. Since, meanwhile, the area under the 
bump remains almost constant, this additional structure 
produces a contribution of order N to the average value of 
Nfi. Consequently, r is finite above the transition. While 
in a realization belonging to the small- A^/^ regime the ru- 
mor remains localized as in the case of p < Pc, sl typical 
realization contributing to the bump includes propaga- 
tion through several shortcuts, thus attaining distant re- 
gions in the system. 
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The localization-propagation transition of our model 
has been described in terms of static features, namely 
the final refractory population Nfi, measured when all 
the interaction events have ceased. In the next section, 
we focus on our central interest here and study the dy- 
namics of the propagation process. 



for an infected individual. Then, at each evolution step, 
either the number of infected individuals n/ increases to 
nj + l at the expense of the susceptible population, or nj 
decreases to n/ — 1 and the refractory population grows 
accordingly. Therefore, the evolution of the number of 
both susceptible and refractory individuals is implicit in 
the evolution of n/. In order to give ni as a function of 
time, it must be taken into account that the number of in- 
fected individuals varies and, consequently, the real-time 
duration of an evolution step changes. At each evolution 
step, in fact, time is to be updated according to 



tn 



(3) 



where the constant to fixes time units. We choose to — 1, 
so that the unit of time can be associated with the typical 
time needed by a single infected individual to establish a 
contact with one of her neighbors. 
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FIG. 1. Normalized frequency distribution /{Nr) of the 
number of refractory individuals at the end of the evolution, 
Nr, for K — 2 and two values of the small- world randomness, 
(a) p — 0.05 and (b) p — 0.3. Different symbols correspond 
to = 10^ (triangles), N = 10* (squares), and iV = lO'' 
(circles). FVequency counts were obtained from series of 10^ 
realizations for each parameter set. 

For p « pc, the distribution (not shown in Fig. ||) ex- 
hibits power-law dependence for moderate values of Nn, 
/{Nr) ^ N^" with a « 1.5. The power-law regime ter- 
minates at a smooth cut-off, whose position shifts to the fi (f^ 
right as the system size increases, approximately as N'^'^ ^ 
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III. EVOLUTION OF THE INFECTED 
POPULATION 

A complete characterization of the propagation process 
in our model is given by the evolution of the infected pop- 
ulation. Initially, all the population is susceptible, except 



FIG. 2. Evolution of the number ni{t) of infected indi- 
viduals as a function of time in two single realizations on a 
10^-node small-world network with K = 2, for (a) p — 0.05 
and (b) p = 0.3. 
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Figure g shows the evolution of the number of infected 
individuals as a function of time for two single realiza- 
tions with = 10^ and K = 2, for two values of the 
network randomness, p — 0.05 and 0.3. For p — 0.05 the 
final number of refractory individuals is Nfj = 32, while 
for p — 0.3 we have Nji — 22, 258. This latter realization 
belongs to the large-A^i? bump structure in /{Nr) for 
the corresponding value of p. Realizations for the same 
randomness but in the small- A^a, region are qualitatively 
similar to that shown in Fig. Ba. 
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FIG. 3. Total time T up to the extinction of the infected 
population, as a function of the final number of refractory 
individuals Nr, for networks with K — 2 and A*' = 10^ (tri- 
angles), 10* (squares), and 10^ (circles). The randomness is 
(a) p = 0.05 and (b) p — 0.3. Both T and Nr were measured 
in 10 realizations for each value of p. In the cases where 
several values of T were obtained for the same value of Nr, 
they were averaged. The straight lines have slope 2/3. The 
insert in (b) shows a close-up of the "cloud" at large values 
of Nr for TV = 10*, in linear scales. 

The graph of nj as a function of time for small p — or, 
for p > Pc, in the region of small Nji — is reminiscent of 
a random walk. As a matter of fact, the evolution of 
the number of infected individuals can be thought of as a 
one-dimensional random walk in the n/-space, starting at 



n/ = 1 and with an absorbing boundary condition at the 
origin, nj = 0, where the evolution terminates. Equiv- 
alently, we may think of a first-passage-time problem, 
with respect to the origin, for a random walker starting 
at 71/ = 1 [Q. This analogy, however, is difficult to ex- 
ploit, since in our case the random walk would be biased 
by a complex time-dependent asymmetry. In fact, the 
probability for nj to grow or decrease depends not only 
on nj itself, but also on the number of both susceptible 
and refractory individuals. The effect of this bias would 
be particularly strong for the large- iV^ realizations with 
p > Pc- In this case, indeed, the evolution of nj{t) does 
not resemble a random trajectory but mimics determin- 
istic dynamics affected by a moderate level of noise (see 
Fig. g,). 

The first-passage-time analogy suggests anyway that 
a compact quantitative characterization of the propaga- 
tion process is given by the total time T elapsed up to the 
extinction of the infected population, and the maximum 
number of infected individuals during the evolution, Nj. 
In the associated random walk, these two quantities cor- 
respond to the first-passage time and the maximum span 
from the origin, respectively. In order to compare with 
our previous results, we measure T and Nj as a func- 
tion of the final refractory population Nr. Note that 
Nr is directly related to the duration of the propaga- 
tion process measured in evolution steps. In fact, since 
the final number of infected individuals is zero, each step 
where a susceptible individual becomes infected must be 
compensated by a step where an infected individual be- 
comes refractory. Since an extra step of this latter kind 
is needed for the first infected individual, the total num- 
ber of steps necessary for the extinction of the infected 
population is exactly 2Nr — 1. On the other hand, due 
to the changing duration of steps in real time, Eq. (||), 
the connection between T and Nr is more complex. 

In Fig. ^ we present measurements of the total time 
T as a function of Nr over series of IC* realizations, for 
SWNs with K ~ 2 and three different sizes, and for two 
values of the randomness p. For small p there is a quite 
well defined power-law dependence, T ^ NJ^, spanning 
almost two orders of magnitude in Nr. Linear fitting of 
these data yields an exponent r close to 2/3. This re- 
sult differs from the value predicted by the random- walk 
analogy for an unbiased random walk in the n/-space, 
which gives r = 1/2. On the other hand, it can satisfac- 
torily reproduced by a random walk with constant bias, 
with probabilities P = 0.6 of moving towards -l-oo and 
1 — P = 0.4 of moving in the opposite direction. We 
recall, however, that the analogy would be strict for a 
time-dependent bias only. For p > pc the above power- 
law dependence is still observed in the small- iV^ regime, 
but apparent deviations appears as Nr grows. In partic- 
ular, the detached "cloud" of dots observed for TV = 10'' 
and 10^ at large values of Nr — which correspond to re- 
alizations in the bump structure observed in /{Nr) (see 
Fig. 1^) — does not satisfy the power-law relation. The to- 
tal evolution times associated with such realizations are 
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considerably below those predicted by an extrapolation 
from the small-iV^ region, and the difference becomes 
larger as the size N grows. Note moreover, from the 
insert in Fig. ||b, that inside the "cloud" there is no ob- 
vious correlation between T and Nji, in contrast with the 
small- iVij regime. These features make it evident that a 
qualitative change in the dynamical behavior occurs be- 
tween the regimes of small and large Nji, as illustrated 
in Fig. |. 
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FIG. 4. Maximum number of infected individuals during 
the evolution, Nj, as a function of the final number of refrac- 
tory individuals, Nn, for networks with K = 2 and A'^ — 10"^ 
(triangles), 10* (squares), and 10^ (circles). The randomness 
is (a) p = 0.05 and (b) p = 0.3. Data were obtained from the 
same numerical realizations as those of Fig. ^. The straight 
lines have slope 1/2. The insert in (b) shows a close-up of the 
"cloud" at large values of Nr for A'^ = 10*, in linear scale. 

Essentially the same features are found for the depen- 
dence on Nji of the maximum number of infected indi- 
viduals during the whole evolution, Nj, shown in Fig. 
^. Now, however, the exponent in the power-law rela- 
tion Nj ~ N^, observed to hold in the small- A^j^ regime, 
is close to 1/2, which does coincide with the result for 
an ordinary unbiased random walk. For large N^, the 
"clouds" of dots quoted above deviate in this case to 



higher values of iV/. A detail of the "cloud" for iV = 10^, 
shown in the insert of Fig. reveals a remnant corre- 
lation between Nj and Njj. 

Note that the power-law dependence of T and Nj on 
Njj, implies the relation T ^ Nj, with ^ = t/i/ w 1/3. 
This relation is expected to hold for small randomness 
or, more generally, for small N^. 

Through the study of the evolution of the infected pop- 
ulation, we have so far examined the dynamical prop- 
erties of our model for just two values of the small- 
world randomness p, below and above the localization- 
propagation transition at Pc- It is now worthwhile to dis- 
cuss how the results change as p is varied. For < p < pc, 
as a matter of fact, the power-law dependence of T and 
Nj on Nn is not modified. The exponents t and v are 
the same within our numerical precision. On the con- 
trary, above the transition, substantial changes affect the 
frequency distribution of the final refractory population 
Njj and its relation to T and Nj. 

In the first place, the relative number of realizations in 
the small- A^jj regime and in the bump at large Nn varies 
considerably with p. Figure |^ shows the fraction p of 
realizations in the large-Nn bump as a function of p for 
three values of TV and K — 2. This fraction grows from 
p ~ 0.35 for p = 0.3 to p fa 0.7 for the maximum random- 
ness p ^ 1. As expected, realizations where the rumor 
attains a significant fraction of the population become 
more frequent as the network randomness grows. Note 
moreover that the dependence with the system size N is 
quite weak, but there is no clear indication of saturation 
for large N. 
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FIG. 5. Fraction p of realizations in the bump at large Nr 
as a function of the small-world randomness p, for networks 
with K ^2 and iV = 10^ (triangles), 10* (squares), and 10^ 
(circles). Data obtained from series of 10^ realizations for 
each values of p. 

As for the dependence of the total evolution time T 
and the maximum number of infected individuals Nj on 
the final refractory population Nji, the small- A^^; regime 



5 



exhibits no significant modifications as p changes. The 
power-law dependence with the same exponents is main- 
tained, as expected from the fact that the reahzations in 
this regime correspond to propagation of the rumor over a 
limited neighborhood of its origin. In contrast, the bump 
in f{Nii) varies in position and the corresponding values 
of T and Nj change. In Fig. ^, wc show the average val- 
ues of T and Nj as a function of the fraction r = N^/N 
corresponding to realizations in the bump for several val- 
ues of p. Roughly speaking, each dot represents the cen- 
ters of the "clouds" referred to in connection with Figs. 
^ and H now for varying randomness. For fixed A^, the 
final refractory population grows with p. For the 
largest systems, in fact, N^^ is practically doubled as p 
varies from 0.3 to 1. Wc thus verify again that propaga- 
tion is more efficient for larger randomness. This effect 
is enhanced by the fact that, at the same time, the maxi- 
mum infected population Nj increases and the total time 
T decreases. For large systems, T is reduced by a factor 
of 2 whereas Nj grows by a factor of 5, approximately. 
In summary, the process becomes simultaneously more 
effective and more rapid. Note, finally, the strong satu- 
ration in the values of r, T, and Nj as the randomness 
approaches its maximum p = 1. 
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FIG. 6. Average values of the total evolution time T 
(empty symbols) and the maximum infected population Ni 
(full symbols) in the large- A'^^ regime as functions of the av- 
erage fraction of final refractory individuals, r = Nr/N, for 
networks with K = 2 and A'^ — 10^ (triangles), 10'' (squares), 
and 10^ (circles). Dashed lines have been drawn as a guide 
to the eye. In each data set, the leftmost and rightmost dots 
correspond to p — 0.3 and p — 1, respectively. From left to 
right, the randomness changes by steps of size 5p — 0.1. For 
A'^ = 10'* and 10"", the value p = 0.35 is also included. Each 
dot stands for an average over 10^ realizations. 

Let us end this section by addressing the effects of 
changing the average number of neighbors per node, 
given by the parameter K. As for the localization- 
propagation transition, a growth in the number of neigh- 



bors implies that the critical randomness Pc decreases 
and that, for fixed p, the final fraction of refractory in- 
dividuals increases [|l^. These two results agree with 
the expected fact that propagation is more efficient for 
larger K. The same trend is observed in the parameters 
that characterize the evolution of the infected popula- 
tion: while the total time T decreases, the maximum 
number of infected individuals Nj grows. Our results for 
K — 2, in any case, are not qualitatively changed when 
other values of K are considered. 



IV. PROPAGATION ON DYNAMIC SMALL 
WORLDS 

Dynamic small worlds (DSWs) have been introduced 
as a variant to SWNs in the frame of a model for activity 
propagation in a system of mobile automata [ p^ . In- 
stead of considering a frozen disordered interaction net- 
work, DSWs admit interactions between any two indi- 
viduals occupying the nodes of a regular lattice. In the 
case of a one-dimensional array with periodic boundary 
conditions, at each interaction event, the partner of an 
individual is chosen with probability 1 — p among its 2K 
nearest neighbors, K clockwise and K counterclockwise. 
With the complementary probability p, the partner is 
chosen at random from the whole lattice. In this way, all 
individuals have the chance to interact with arbitrarily 
distant partners, but the probability of distant interac- 
tions is controlled by the "randomness" p. 

The change from SWNs to DSWs, which conveys the 
replacement of frozen disorder by a stochastic process, 
is qualitatively similar to the introduction of the so- 
called annealed approximation in the study of disordered 
Boolean (Kauffman) networks A considerable ad- 

vantage of DSWs over SWNs regards the numerical im- 
plementation, which does not require the generation of a 
new lattice at each realization — a highly time-consuming 
step in our specific system. However, the main virtue of 
DSWs — shared with the annealed model for Kauffman 
networks — is that, in principle, they admit a simpler an- 
alytical treatment. In particular, the limit p = 1 should 
be exactly described, in asymptotically large systems, by 
a mean-field-like approach. 

We show in the following that the behavior of the ru- 
mor propagation model on a DSW bears remarkable sim- 
ilarity with the case of a SWN, though some significant 
quantitative differences are detected. Let us first of all 
point out that, as for the dependence of our model on the 
system size A'^ and on the average number of neighbors 
per individual, 2K, the features described in the previous 
section for SWNs are qualitatively reproduced in DSWs. 
Consequently, we do not repeat the analysis for varying 
N and K, and focus here on the specific case A^ = 10^, 
K ~ 2. The behavior for other values of A^ and K can 
be inferred from this case and the results for SWNs. 

It has already been advanced ||l^ that our model 
on a DSW undergoes the same kind of localization- 
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propagation transition found on SWNs. This is illus- 
trated in Fig. ^ where we show the final fraction of 
refractory individuals r — Nn/N as a function of p, for 
N = 10'* and K = 2. The critical point has considerably 
decreased, to Pc ~ 0.06. Meanwhile, as expected, the 
fraction r approaches the solution to Eq. (|l|), r* « 0.796 
as p ^ 1. Note in fact that the original version of the 
model, discussed in Sect. |l[ is a kind of mean-field ap- 
proximation of the small-world case, which becomes ex- 
act for p ~ I. 
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FIG. 7. Final fraction of refractory individuals, 

r = Nr/N, as a function of the "randomness" p on a dy- 
namic small world with total population = 10* and K = 2. 
Each dot stands for an average over 10** realizations. The 
dashed line is a spline approximation, drawn as a guide to 
the eye. 

As for SWNs, the nature of the localization- 
propagation transition in DSWs is revealed by the fre- 
quency distribution of the final number of refractory in- 
dividuals, f{N]i). Figure ^ shows this distribution for 
three values of p. For p = 0.02 we find a rapidly de- 
caying function which, as in the subcritical regime on 
SWNs, results to be roughly exponential. In this case, 
in fact, the contribution of distant interactions is negligi- 
ble, so that no significant differences are to be expected 
between DSWs and SWNs. For p — 0.06, which approxi- 
mately corresponds to the critical point Pc in DSWs, the 
distribution is a power law over a substantial interval, 
/{Nr) ^ A^/j"- Remarkably, the exponent of this power 
law coincides — up to the numerical precision — with that 
obtained at the critical randomness in SWNs, a « 1.5 (cf. 
Sect. Finally, above the critical point {p = 0.1), we 
find that the by-now familiar bump structure at large Nr 
has developed. The strong similarity with the scenario 
on SWNs convincingly suggests that the origin of the 
localization-propagation transition is the same for both 
systems. 

The insert of Fig. || shows the fraction p of realizations 
that contribute to the bump, as a function of p. Compar- 



ing with Fig. ^, which shows the same results for SWNs, 
we note that — apart from the obvious consequences of the 
shift of Pc to the left — the fraction p attains considerably 
larger values. In particular, we find p w 1 for 0.6 < p. 
For such values of p, therefore, the rumor propagates to 
distant regions and attains a finite portion of the system 
in virtually all realizations. An important contribution 
to this difference with SWNs is given by the following 
fact. In our SWNs, links between individuals are bidi- 
rectional. This implies that if an infected individual i 
transmits the rumor to a susceptible neighbor j, there 
is a relatively high probability that, in the future, the 
now infected individual j will (unsuccessfully) attempt 
to transmit the rumor to i and will become refractory. 
These unsuccessful trials for backward propagation are 
by far more improbable in DSWs, especially for large p, 
and the rumor spreading is consequently enhanced. 




FIG. 8. Normalized frequency distribution /(A^a) of the 
final number of refractory individuals Nr on a dynamic small 
world with a total population A'^ = 10* and K = 2, for 
p = 0.02 (empty dots), p = 0.06 (crosses), and p = 0.1 (full 
dots). The frequency distribution was obtained from series of 
10® realizations for each "randomness." The dashed straight 
line has slope —1.5. The insert shows the fraction p of real- 
izations belonging to the large- A'^h bump as a function of p, 
averaged over 10^ reahzations. 

For p < Pc the total evolution time T and the maxi- 
mum number of infected individuals Nj satisfy the same 
power-law dependence of Nr as in SWNs. In fact, as far 
as long-range interactions remain infrequent, the evolu- 
tion on SWNs and DSWs is essentially equivalent. The 
same argument can be extended for p > pc m small- A'^h 
realizations. In this case, however, the exponent r in the 
relation T ^ NJ^ results to be smaller than for SWNs; 
now, we find r « 0.57. The moderate contribution of 
distant contacts is here enough to produce a consider- 
able decrease of the total evolution time. Since contacts 
between any two individuals are now possible, propaga- 
tion on DSWs is faster. 
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FIG. 9. (a) The total evolution time T and (b) the maxi- 
mum number of infected individuals Nj, as functions of the 
final refractory population Nr, on a DSW with N = 10* and 
K = 2, for p = 0.5 (empty dots), p = 0.7 (crosses), and 
p = 0.9 (full dots). Only the large- A''ij region, corresponding 
to the bump structure, is shown. Data were extracted from 
10^ realizations for each value of p. In the cases where sev- 
eral values of T and Ni were obtained for the same value of 
Nr, they were averaged. Square symbols stand for the values 
of Ni and Nr, and the lower bound for T predicted by the 
mean-field approximation. The curves correspond to these 
same values for varying A^. 

The fact that propagation on DSWs is faster and more 
effective than in SWNs becomes dramatically emphasized 
as soon as the lavge-Nn regime is analyzed. Figure ^ 
shows T and Nj as functions of Nji for large "random- 
ness," ranging from p = 0.5 to 0.9. These plots show the 
"clouds" corresponding to the large- iVj^ bump structure, 
and are therefore analogous to the inserts of Figs. ||b and 
|b. The typical values of Nr, T, and Nj in these DSW 
realizations are to be compared with the corresponding 
values for SWNs, shown in Fig. ^. While, for N — 10'', 
the fraction r = Nr/N attains in SWNs a maximum av- 
erage level of about 0.44 for p = 1, in DSWs r reaches 
a typical value close to 0.8 for p = 0.9. The differences 
in T and Nj are even more drastic. In SWNs, their ex- 
treme average values are T « 75 and Nj 440, whereas 



in DSWs they change to T « 25 and Ni ?s 3000. 

As the "randomness" of DSWs approaches its maxi- 
mum value p = 1, the system should be satisfactorily 
described by a mean-field approximation. The analyt- 
ical treatment of our model within such approximation 
is developed in the appendix. There we show that, for 
a given size N, it is possible to predict the values of 
Nj, Nfj, and a lower bound for T. These values are 
compared with numerical results in Fig. ^. The agree- 
ment with the average of numerical results is very good 
for the largest "randomness," p — 0.9. These results, 
however, show considerable fluctuations with respect to 
the values predicted by the mean-field approximation. 
Moreover, fluctuations in Nr, T, and Nj are closely cor- 
related. Note that the same kind of correlations were 
suggested in SWNs by the results shown in the inserts 
of Figs. ^ and Since, as discussed above, DSWs 
can be more efficiently implemented in numerical exper- 
iments, our results in Fig. ^ correspond to a number of 
realizations considerably larger than for SWNs, and such 
correlations become apparent. 

It turns out that the correlations between the values of 
T , Nj, and Nji obtained in single realizations can be ex- 
plained in terms of the mean-field approximation. In fact, 
calculating Nj and Nji, and the lower bound for T from 
this approximation for different values of the system size 
TV— ranging from N « 9200 to 10300— we obtain the val- 
ues shown in Fig. ^ as curves. These values successfully 
reproduce the correlation between the three quantities. 
Specifically, as N increases, T decreases slowly and Nji 
grows, while — as shown in the appendix — Nj and Nr are 
linearly correlated. From a phenomenological viewpoint, 
these results can be interpreted as if in each individual re- 
alization the system appears to have an "effective" size — 
close to, but different from, its actual size — plausibly de- 
termined by variations in the effectiveness with which the 
rumor spreads over the population. 



V. CONCLUSION 

We have here studied the evolution of an epidemic-like 
model evolving on small-world geometries. The dynam- 
ics, which can be interpreted as the spreading of a ru- 
mor, is known to exhibit a transition between regimes 
of localization and propagation at a finite randomness of 
the underlying disordered geometry |10|. Epidemiologi- 
cal models on geometries that plausibly represent social 
networks and information webs — such as small- world and 
scale-free networks [|7| — have recently attracted much 
attention, in view of their potential role in the descrip- 
tion of actual risk situations associated with infectious 
diseases and computer viruses [p|,p^-^. 

In our model, the effectiveness of propagation is char- 
acterized by the total number Nr^ of individuals that have 
been infected during the whole evolution. Generally, in 
a single realization of the process on an asymptotically 
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large system of size N, propagation affects either a van- 
ishingly small fraction of the population, Nn/N 0, or 
a finite fraction r. For a randomness p below the critical 
point Pc of the localization-propagation transition, only 
the small- A''/} regime is observed, and the rumor remains 
localized within a limited neighborhood of its origin. For 
p > Pc, on the other hand, realizations in both regimes 
are observed. The fraction of realizations in the large- A'^r 
regime, in fact, grows as the randomness increases. 

The dynamics of our model is completely described by 
the evolution of the number ri/ of infected individuals. 
A compact characterization of this evolution is given by 
the total time T elapsed up to the extinction of the in- 
fected population, the maximum number Nj of infected 
individuals at a given time, and the total number Nr 
of infected individuals during the whole evolution. The 
effectiveness of propagation increases when T decreases, 
because spreading is faster, and when Nr and Nj grow, 
because the rumor reaches a larger population. 

Our results for small-world networks can be summa- 
rized as follows. For any value of p, the small- A'fl regime 
is characterized by power-law correlations between T, 
Nji, and Nj. It has been suggested that these correlation 
could be explained in terms of a random-walk picture of 
the propagation process in the nj-space. A rigorous anal- 
ogy, however, can only be achieved in terms of a biased 
random walk with a rather complicated time-dependent 
bias. In the large-A^_R regime, the values of T, Nr, and 
Nj obtained in single realizations are distributed around 
certain typical values, which vary as the network ran- 
domness p changes. Specifically, as p grows, T decreases 
and both Nr and Nj increase, indicating that the propa- 
gation process becomes increasingly effective. The three 
quantities show a quite marked saturation as the random- 
ness approaches its limiting value p = I. The effective- 
ness of propagation is also improved, as expected, when 
the average number of neighbors per individual grows. 

We have also studied these features in a so-called dy- 
namic small world. Instead of considering a frozen net- 
work of interaction links, dynamics small worlds admit 
that distant contacts can occur between any two individ- 
uals, chosen at random at each evolution step. We have 
here shown that, as far as our model is concerned, propa- 
gation in a dynamic small world is qualitatively the same 
as on a small-world network. Namely, the same kind of 
correlations between T, Nr, and Nj and the same de- 
pendence with p observed on small-world networks are 
reproduced in dynamic small worlds. The main quanti- 
tative difference between both cases is that in dynamics 
small worlds the effectiveness of propagation is consider- 
ably higher. Evolution times are overall shorter and in- 
fected populations larger than on small-world networks. 
In qualitative terms, this is plausibly due to the fact that, 
in dynamics small worlds, the average number of interac- 
tion partners per individual is very large and the effect of 
backpropagation is comparatively negligible, especially, 
in large populations. 

To our knowledge, this is the first time that the evo- 



lution of a dynamical process is compared in detail on 
small-world networks and in dynamic small worlds. It 
may be conjectured that our main conclusions regarding 
this comparison will hold for a large class of processes. A 
systematic comparison would in fact be desirable since, 
though small- world networks have attracted considerably 
more attention than dynamic small worlds, the latter 
have the advantage of easier analytical and numerical 
treatment and, moreover, provide a more realistic model 
of social systems with mobile individuals. 
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APPENDIX: MEAN-FIELD APPROXIMATION 

Taking into account that the susceptible, infected, and 
refractory populations — which we denote ns{t), ni{t), 
andnij(t), respectively — satisiy ns{t)+ni{t)+nR{t) = N 
for all times, the mean-field evolution of our model is 
given by the two equations 



and 



hi — nj \ 1 



ur = ni- 



' N 



ni + ur 



(Al) 



(A2) 



These equations can be implicitly solved in terms of the 
auxiliary variable 



ni{t')dt'. 



(A3) 



We point out that the introduction of the variable s in 
this continuous approximation is fully equivalent to the 
change from the real time scale to the measure of time 
in evolution steps used in the discrete model. 

With the initial conditions nj(Q) — 1 and nR{Q) — 0, 
the solutions to Eqs. (Al) and (|A2|) are 



ni{s) = 1 - s + 2(iV - 1)[1 - exp(-s/7V)] (A4) 



and 



nR{s) = s - (iV - 1)[1 - exp(-s/Ar)]. 



(A5) 



From these solutions, some relevant quantities can be 
immediately calculated. The total number of steps S 
needed for the extinction of the infected population is 
the positive solution to nj = 0, which corresponds to the 
transcendental equation 



S -1 = 2{N exp(-5/iV)]. 



(A6) 



For each value of N, this equation can be accurately 
solved by numerical means. For asymptotically large N, 



it can be shown that S = kN, where k 1.594 is the pos- 
itive solution to k — 2[1 — exp(— fc)] [note that k = 2r*; 
cf. Eq. (|l])]. The final number of refr actory individuals, 
Nr, can be evaluated from Eq. (A5) for s = S. This 
yields 



Nr = 



S+1 



(A7) 



Actually, this result holds not only in the mean-field ap- 
proximation, but for any value of p and N in both SWNs 
and DSWs, as discussed in the main text. 

The step sj at which the infected population nj(t) 
attains its maximum Nj is given by hi = 0, i. e. 



si = iVln[2(iV- l)/iV]. Replacing in Eq. ( |A4[ ) we obtain 



Nr = N ■ 



Nln 



2{N - 1) 
N 



(l-ln2)7V, (A8) 



where the right-hand side approximation holds for large 
N. In this limit, combination of the above results makes 
it possible to show that 



Nr 



2(1 - In 2) 



-Nr « 0.385iVH. 



(A9) 



The only problematic point in the comparison of 
the mean-field approximation with the original discrete 
model is the evaluation of the total (real) time T elapsed 
up to the extinction of the infected population. As a 
matter of fact, in the mean- field approximation nj(t) 
decreases asymptotically and vanishes only for t — > oo 
(note that, however, this limit corresponds to a finite to- 
tal number of steps S). The total time T must therefore 
result from a plausible definition using the mean-field re- 
sults. In the main text we use the following criterion. We 
define T as the time needed for nj(t) to attain again its 
original value nj = 1. In other words, T is the nontrivial 
solution to ni(t) = 1, which can be acc uratel y ob tained 
from numerical integration of Eqs. (Al) and (|A2|). This 
gives a lower bound for the actual total evolution time, 
that can be directly compared with our numerical results. 
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